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Abstract 

We shall investigate the possibility of formulation of varying speed of 
light (VSL) in the framework of Palatini non-linear Ricci scalar and Ricci 
squared theories. Different speeds of light including the causal structure 
constant, electromagnetic, and gravitational wave speeds are discussed. 
We shall see that two local frames are distinguishable and discuss about 
the velocity of light in these two frames. We shall investigate which one 
of these local frames is inertial. 



1 Introduction 

There are many theoretical and experimental suggestions that the values of 
fundamental physical constants are not actually constant and may vary in space 
and time pQ. One of these constants which has attracted considerable attention 
recently [2] is the speed of light. The first modern VSL theory was Moffat's 
ground breaking paper [3], where spontaneous breaking of Lorentz symmetry 
leads to VSL theories. Then these theories followed by Albrecht, Magueijo and 
Barrow [4l [5] as an alternative to the inflation mechanism to solve some puzzles 
of Big-Bang cosmological models[H |B]. In such a VSL cosmology the Lorentz 
invariance is explicitly broken and there is a preferred reference frame in which 
the physical laws should be formulated. They called this frame the cosmological 
frame and assumed that in this frame interactions have to be introduced via a 
minimal coupling principle according to which the curvature tensor is computed 
ignoring the variation of the speed of light. In this way the field equations are 
Einstein's equations with the varying speed of light appearing in the right hand 
side. 

Generally in VSL models either a pre-set function for the speed of light is 
considered [5, 7 J or there is a term in the Lagrangian determining the dynamics 
of the speed of light [8]. However Magueijo in [8] has proposed a generalized 
varying speed of light theory preserving the general covariance and local Lorentz 
invariance. This is achieved by introducing a time-like coordinate x° which is 
not necessarily equal to ct, and the physical time t can only be defined when 
dx° / c is integrable. 
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It has to be emphasized here that there are different constants which can be 
interpreted as the velocity of light. As Ellis and Uzan[9j have shown one have to 
distinguish between cem (the electromagnetic wave velocity), cst (the space- 
time causal structure constant), cqw (the gravitational wave velocity) and ce 
(the space-time-matter coupling constant appearing at the right hand side of 
Einstein's equations). They have shown that assuming the standard Lagrangian 
of the electromagnetism and general relativity, one has cem = cst and also that 
c gw = cst- They have also shown that having the correct Newtonian limit leads 
to c E = cst- 

As a result, in dealing with any modified gravity and/or electromagnetism 
theory one should re-examine the meaning and the relation between these con- 
cepts of the speed of light. Also to make the concept of varying speed of light 
meaningful, one should consider dimensionless quantities like ce/cq, cqw/cq, 
cem I cq, and cst/cq, where cq is a constant of dimension velocity and we can 
choose it equal to 3 x 10 8 meters/ seconds in the MKS units. In this way for the 
standard theory of electromagnetism and gravity we have: ce/cq = cgw/cq = 
cem/ cq = c S t/ cq = 1. 

On the other hand, there are many different motivations to replace the stan- 
dard general relativity by a more complete one which has richer space-time con- 
cepts. The existence of solutions with curvature singularities or closed-time-like 
loops in the standard general relativity are some of these reasons. Furthermore 
the flat rotation curves of galaxies and the current accelerated expansion of the 
universe introduce dark matter and dark energy with an amount about 96% 
of the contents of the universe may be interpreted as signals for need for an 
extended theory. 

There are different approaches to modify Einstein's gravity theory. These 
can be put into different classes [10] . Scalar-Tensor theories [TT], higher di- 
mensions theories (Brane- World scenarios[12] and Kaluza-Klein theory [13]) and 
theories with a modified Lagrangian (f(TZ, 1Z a b, • ■ ■)) in both Palatini and metric- 
affine formulations [TH Q~5] , are some of them. 

In the Palatini gravity, the metric and the connection are considered as 
different degrees of freedom, but the matter Lagrangian does not depends on 
the connection. For the case of the standard gravitational Lagrangian, y/~glZ, 
this theory is identical with the Einstein's theory. If the matter Lagrangian 
depends on the connection then the theory is called metric-affine gravity. Two 
special cases in which the Lagrangian is y/—gf(TZ) (non- linear Ricci scalar) and 
y/^glZ + \/—gf(TZ ab TZ a b) (non-linear Ricci squared) are considered frequently 
in the literature and have some interesting results as well as some drawbacks 

01 Eg. 

Since in such theories the connection is different from the Christoffel sym- 
bols, one can introduce two different local frames. One in which the connection 
is locally zero and the other in which the Christoffel symbols are locally van- 
ishing. These two frames are conformally related for Palatini f(TZ) gravity and 
thus the causal structure of both frames is the same, but this does not mean 
that they are completely equivalent from the physical point of view. We shall 
discuss about this point later. For Palatini f(lZ ab lZ a b) gravity, these are not 
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conformally related and thus the causal structure of these two frames is differ- 
ent. Investigation of the speed of light in these two frames helps one to choose 
one of these local frames as the local inertial frame. 

In this paper after a brief review of the Palatini f(lZ) and f(lZ ab lZ a b) gravity, 
we shall discuss about two possible local frames and about choosing one as the 
local inertial frame. Then, we shall evaluate the different velocities of light, 
including the gravitational and electromagnetic wave velocities, and the space- 
time causal structure constant. 

It has to be noted here that according to Ostrogradski theorem[T7], higher 
derivative Lagrangians may suffer from instability problems. According to Os- 
trogradski, Hamiltonians associated with Lagrangians containing more than first 
time derivative would have linear instabilities, provided that higher derivatives 
cannot be eliminated by partial integration. It can be shown such instabilities 
does not exist for f(lZ) metric theories [18]. 

In the Palatini formulation of f(lZ) and f(lZ ab lZ a b) theories, the Lagrangian 
depends only upon the first derivative of the fields and thus the result of Ostro- 
gradski theorem is not applicable here. 

This is the first reason why the metric formulation is not used here. There 
is also another reason. In the metric formulation, since the connection is equal 
to the Christoffel symbols there is only one local frame. This, as it would be 
clear at the end of this paper, leads to a constant speed of light. 

2 Non- linear Ricci scalar and Ricci squared Pala- 
tini gravity in a nutshell 

The astronomical evidence for an accelerating cosmic expansion has stimulated 
many investigations into the other possible deviant gravitational effects, which 
might be responsible for this unexpected dynamics[l9j. Many investigators have 
also focused their attentions on modifying general relativity in the large scales. 
One example of them is provided by f(TZ) gravity models and the other example 
is f(lZ,lZ ab ) gravity models. 

When the Einstein-Hilbert action is modified by adding some general func- 
tion of 1Z or !Z ab lZ a b, it becomes necessary to distinguish between two different 
variational approaches for deriving the field equations. In the metric approach, 
as in Refs IS] , the metric components g a b are the only dynamical quanti- 
ties and the field equations are generally of fourth-order. Within the Palatini 
variational approach, on the other hand, one assumes that the space-time is 
described by two independent geometrical objects, the metric and the connec- 
tion. The first one defines the causal structure of the space-time and the latter 
defines the affine structure of the manifold. As a result the equations of motion 
are of the second order. 
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2.1 Palatini non-linear Ricci scalar gravity 

Let us first start with the non-linear Ricci scalar gravity in the Palatini formal- 
ism. The Lagrangian density of the theory is chosen to be an arbitrary function 
of the scalar curvature, f(TZ). According to the above discussions, the action of 
the Palatini f(TZ) theory is: 



(1) 



where K = 8itG/cq, and A m is the matter action. Note that in the Palatini 
formalism the matter action does not depend on the connection. Since the 
affine connection is different from the Christoffel symbols, we can define two 
kinds of derivatives as: 



V a X b = d a X l 



x< 



(2) 



D a x b = d a x b + r b ca x 



(3) 

Variation of the action ((!]) with respect to the metric (g a b) and the connection 



(r£ c ) leads to: 



SA = 



2 k 



d i X\/—g 



f'K ab - y 9ab - nT ab ) Sg ab 



f'g ab (D c ST c ab -D b ST: c )} 



(4) 



in which a prime denotes differentiation with respect to 1Z and T ab is the matter 
energy-momentum tensor defined as: 



6 A m =~\j d 4 xV=gT ab 5g ab 



(5) 



Considering the metric and the connection as independent variables one arrives 
at the following equations of motion: 



f'(TZ)TZ ab ~ ^f(TZ)g ab = nT ab 



and 



d c {V=gfg ab ) = o 

The last equation shows that making a conformal transformation as: 

Kb = f'gab (8) 
leads to the fact that the connection is compatible with the metric h ab given by: 

1 



(6) 
(7) 



lab 



2f> 



2S c (a d b) f - g ab g cd d d f 



(9) 
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A very interesting result can be obtained by investigation of the general co- 
variance of the action (fTJ) . Let us make an infinitesimal transformation of the 
coordinates x a — > x a + e a . Using the relation ([5]) and the fact that under coordi- 
nate transformation we have Sg ab = g< a d c e b \ after some ordinary calculations, 
one finds that the corresponding Noether equation is: 

V Q T Qb = (10) 

For an explicit derivation of this relation the reader is referred to [22j . The above 
equation means that the covariant divergence of the energy-momentum tensor 
does not vanish. This is not a bad signal since one can argue that it should be so. 
Any conservation relation should be of the form V a j a = —j=d a (\/~gj a ) = 0, 
to lead to the conservation of J d 3 Xy^—gj a n a (with n a a time-like vector). Note 
that the correct measure for this integral is y/—g and not \f—h, since this is 
g a b which is used to measure distance, area and volume, as we are dealing with 
a space-time structure in which the metric is g a b and the connection is T c ab . 
Extension to the case of higher rank tensors is straightforward. 



2.2 Palatini non-linear Ricci squared gravity 



In order to make a more general extension to the standard gravity one may add 
to the Einstein-Hilbert action a function of both the Ricci scalar and the Ricci 
tensor. A simple case that we shall use here is that this extra term is only a 
function of Ricci tensor, i.e.: 



A 



d^Xy 



K + f(K ab K ab 



2 k 



(11) 



By varying the action (jlip we have 
1 



8A 



d 4 x- 



2k 



g K ab + 2FK c a K bc - -g ab [K + f(K ab K ab )] - KT ab 5g 



ab 



d A x^^-g{g ab + 2FTZ ab )[D c 6T c ab - D b ST c ac 
Zk 



(12) 



where F = df/dS with S = 1Z lZ ab - Therefore the modified Einstein equations 
is: 

J „ , ^ (13) 



TZ ab + 2FTZ c a TZ bc - -g ab [n + /] = KT ab 

On the other hand, variation of the action with respect to the r^ c leads to: 

D e [V=9(g ab + 2Fg ac K cd g bd )] = (14) 

Just like in the Palatini f(lZ) models, this equation implies some relation be- 
tween the physical metric g ab and a new metric (G ab ) whose Christoffel symbol 
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Since finding the solutions of the eqs. p^]) , (fl"4"|) and then finding the relation 
between Gab and g a b is generally a difficult task, one can use the covariant 1 + 3 
approach (see [25J[MJ[55]) at least for the cosmological solutions. 

Let's to assume that there is a family of preferred worldlines representing 
the average motion of matter at each point (as is the case for cosmology). This 
is the 4— velocity u a = 2§— normalized as u a u a = 1. We assume this 4— velocity 
is unique: that is, there is a well-defined preferred motion of matter at each 
space-time event. Given this, two projection tensors can be defined. The first 
is U£ = u a Ub which projects anything parallel to the 4— velocity vector, and the 
second is h a b = gab — UaUb, which determines the (orthogonal) metric properties 
of the instantaneous rest-spaces of observers moving with 4— velocity u a . 

The energy-momentum tensor T a b can be decomposed relative to u a in the 
form 

Tab = pu a u b + 2q( a u b ) - phab + n ab (15) 

in which , ir a b — h^h^Td+phab is the projected symmetric trace free anisotropic 
pressure (stress), p — T a bU a u b is the relativistic energy density relative to u a , 
q a = hfu c T c d is the relativistic momentum density and p = T a b is the 

isotropic pressure. 

In a similar way lZ a b can be written as: 

TZab = Aw a ii 6 + Eh a b + 2w (a T b) + E ab (16) 

Substituting cqs. flT)]) , (HHJ) into cq. (fT3"|) , we get the following four equations, that 
determines the coefficients 

A + 2FA 2 -i(A + 3S + /) = K p (17) 

S + 2FS 2 -i(A + 3S + /) = -«p (18) 
[l + 2F(A + S)]r a = K g Q (19) 

(1 + 4FH)S o6 = K7T ab (20) 

where f, F are functions of TZ ah TZ a b = A 2 + 3S 2 .Thus given the form of / and 
the values of p,p, q a , i^ab the quantities A, H, T a , ir a b, can be obtained from the 
above equations at least numerically. So we have TZ ab from eq. (|16|) . 

Now we should solve the eq. flT4l . As it is stated previously one can introduce 
a new metric Q a b by means of the following relation 

bd\ 



V^GG ab = V=g(g ab + 2F 9 ac n cd g bd ) (21) 

After some calculations we can find the relations between the two metrics Q ab 
and g a b and their inverses as 

Gab = ^9ab + Cab, (22) 

G ab = \g ab + C h (23) 
A 
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where 



A = V(1 + 2 J FA)(1 + 2FS) (24) 
1 + 2FS 



1 + 2FA 



(25) 



2F 

t ab = X(uj - l)u a u b - A^/ujFu (a T b) j=T, ab (26) 

C ab = v ( - - A u a u b + -^^= [2 J a T b ) + S ah ] (27) 
The relation between two different Christoffel symbol is 

a -{£}+* <*> 

where 

7 h a c = ^KVcA + <5 c Q V b A - g bc V a \] + ^[C 6 °V C A + C c Q V b A - g bc Q ad V d \] 

+ ^bC. + v c 4 a - v Q £ 6c ] + ^C ad [v c & d + v fc e C d - V d 6c] (29) 

and also we have 

TZ ab = R ab + V cl C ab - ^bY ac + lib-fed - laclbd (30) 

in which R ab is the Ricci tensor with respect to Christoffel symbols. For more 
details see [26l [16] . 

For our later use, we can write the modified Einstein eq. (|13p as 

Rab - TjgabR = KT ab + nT^ f (31) 



where 
and 



K T e J f = ^g ab (f + g cd U cd ) - n ab - 2FK a K cb (32) 



n Qb = v cl c ab - x/ bl c ac + _ (33) 

As a special case for the FRW cosmological models we have q a = Tr ab = 
so that T ab — pu a u b — ph ab . According to this and eqs. (fl"9|) . (|20[) we have 
T a = S ab = 0. The remaining parameters satisfy the following equations: 

A + 2FA 2 -i(A + 3S + /) = K p (34) 
S + 2FS 2 - -(A + 35 + /) = -kp (35) 

and 

U = \(u - l)u a u b , C b = U- - l)u a u b (36) 

A to 
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So we have 



Q' 




(37) 



(38) 



It should be noted at this end that the conservation relation of the energy- 
momentum tensor is again of the form V Q T ah = 0. This is because of the fact 
that the gravitational and matter parts of the action functional are scalars under 
general coordinate transformations. And as it is proved previously this leads to 
the above conservation law. 

3 Local inertial frame 

In order to discuss some kinds of the speed of light, we need to go to the 
local frame. For example for defining cst one should examine the local causal 
structure. Also to separate the gravitational effects from the effects caused by 
variation of the speed of light in the cem it is helpful to go to the local inertial 
frame. 

In the Einstein's theory of gravity one has the equivalence principle which 
asserts that there is some local frame in which gravity is locally absent. This 
is hard-coded in the theory such that the local frame in which the metric is 
Minkowskian has zero connection (the Christoffcl symbols). In the Palatini 
formalism as while as f(TZ) — 1Z, the same is true, but for a general f(1Z) this 
fails. That is to say, one is able to define two local frames. In the first one, 
metric is locally Minkowskian, but the connection is not zero. In fact we have: 
9ab = Va,b, {ab\ = ^ anc ^ ^ab = 7afc- There is also a second local frame in which 
the connection is zero and thus: g a b — jrVab, {^ b } = —lab an d F ab = 0. In 
order to understand the physical significance of the above mentioned frames, 
let's to investigate the trajectory of a test particle in these two frames. To do 
so, we consider a dust with T ab — pu a u b , where p is particle density and u a is 
the velocity four-vector. The conservation relation (JTUJ) leads to: 



Note that according to the discussion of the previous section, this conservation 
equation should be evaluated by V. Combination of these two equations leads 
to the trajectory of a test particle which is u a V a u = 0, or: 



Va(pu a U b ) = 



(39) 



and we have to consider the particle conservation law, too: 



V„(pu a ) = 



(40) 




(41) 



In the first local frame, the particle trajectory is a straight line 



while in the second frame one has 



d 2 x' 



dx a din/' 



+ ^V ab d b f = 



(43) 



ds 2 



ds ds 



Note that the same is true for non-linear Ricci squared theories, except that 



the form of jg c is different. Therefore, generally the trajectory of a test particle 
in the second frame is 



Now the question is that which one of these frames is the local inertial 
frame? It has to be noted that although in the first frame the particle moves on 
a straight line but this does not mean that equivalence principle is satisfied since 
gravity is present via the non-vanishing connection in other physical relations. 
Although we are talking about Palatini formalism in which the matter action is 
not coupled to the connection, but this is not the whole story. 

In Palatini f (Jl ab 1l ab ) and f(TZ) gravity we assumed that the space-time 
structure is defined via the metric (as the device for defining length) and the 
connection (as the device for parallel transporting). This means that in any 
experiments that the concept of parallelism is important the effect of the con- 
nection can be seen. 

For example suppose that we want to measure the focal length of a lens in 
the first local frame. To do so we need to send at least two rays of light towards 
the lens. Although in this frame the Maxwell equations has no dependence on 
the connection and thus the light rays behave like in special relativity, but to 
make them parallel we need to use the connection which is not zero. In other 
words in the first local frame gravity is present in defining the focal length of 
the lens. The distance between geodesies and the description of a congruence 
of them are of course given by the affine connection. 

To state this important point in another way, let us to stress that the theory 
is not only the field equations. It is the field equations derived from the action 
defined on some spacetime with predefined properties. For the Palatini theory 
one assumes the spacetime has an independent affine connection and thus any 
parallel transport should be evaluated using it. So one expects to have changes 
in the geometrical concepts like geodesic deviation and Raychaudhuri's equation 
representing how a flux of geodesies expands. In this theory the existence of two 
different connection fields has new consequences. The geodesies are determined 
by the Christoffel symbols (this choice is motivated by energy-momentum con- 
servation) but the equation that governs the evolution of the deviation vector 
involves the affine connection (motivated by the fact that the covariant deriva- 
tive or parallel propagation along any arbitrary curve is defined by the affine 
connection). For more details see [27] . 

Therefore it seems that the first local frame is not the local inertial frame. 
On the other hand, in the second frame the space-time for Ricci scalar model is 
locally conformally flat and the connection is zero and for Ricci squared model 
only the connection is zero. Since the Christoffel symbols arc not zero, mat- 
ter fields and the test particle trajectory are not the special relativistic ones, 




(44) 
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it is given by equation In order to understand physical implications of 

this equation, it is instructive to define another concept of the speed of light, 
the speed used to synchronize the clocks. In the standard relativity, the syn- 
chronization procedure is as follows. In the local inertial frame, point 1 sends 
an information signal (with velocity cc — Co) to point 2 at a distance d£ from 
point 1. Then if the observer at point 2 adjust its clock d£/co ahead, the two 
points are synchronized. As a result the proper time in this frame is defined as 
ds = ccdr. It is now instructive to write the test particle trajectory (equation 
(@1])) in terms of the proper time (dr = ds/cc)- 

J 2 -r a , dx b dx c dx a 1 dc c , , 

Ibc-rr-rr - -tt — -rr =0 ( 45 ) 



dr 2 c dr dr dr cc dr 

Therefore, although in this frame the connection is zero and thus can be 
interpreted as the local inertial frame, but the particle's trajectory is not a line. 
One can identify the extra terms with the force due to variation of the speed of 
light. To see this, consider the case of Palatini f(R) gravity, in the second local 
frame. The observer at point 2 has to adjust it clock ahead with the amount 
\fJ 1 dljcQ (with d£ 2 = f'\dx\ 2 ). This show that for synchronization one should 
use the velocity cc = co/v 7 / 7 - Therefore the test particle motion is given by: 

dx a 1 dc C ab 2®b C C 



dr 2 dr cc dr cc 

It is clear now that the extra terms are forces exerted on the particle because 
of the variation of the clock synchronization speed of light. A similar discussion 
can be done for non-linear Ricci squared theories. 



4 Speed of Light 

Let us now consider the meaning and the properties of the speed of light in 
the framework of the geometrical structure of the Palatini f(TZ) and f(TZ ab TZ a b) 
gravity theories. As it is discussed earlier, one should distinguish between dif- 
ferent velocities of light. 

4.1 The space— time— matter coupling constant 

As it is clear from the actions chosen for both nonlinear Ricci scalar and Ricci 
squared Palatini theories, we have assumed that ce = cq. Any attempt to 
make this velocity varying should be done by making it dynamical and adding 
dynamical terms for it to the action. This leads to some scalar-tensor theory 
and is discussed in the literature [8l l9l [28] . 

4.2 The gravitational wave velocity 

In order to drive the gravitational wave velocity, one needs to linearize the field 
equations around the no-gravity solution. Since it is not clear what is and 
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whether the linearized solution exists for a general f(lZ) Palatini theorv [T4"l,l23l 
[30] . we restrict ourself to the case in which f(TZ) has not negative powers of the 
Ricci scalar. That is we assume f(TZ) =1Z + a{R 2 + a 2 lZ 3 + ■ ■ ■ where a^s are 
constants. For f(lZ) theory, taking the trace of equation©, we get 

f'(1l)TZ-2f(TZ) = KT (47) 

Assuming T = 0, which includes vacuum and electro-vacuum, for a given /, this 
is an algebraic equation for 1Z. The Ricci scalar will therefore be a constant and 
a root of the equation 

f\Tl )n - 2/(72o) = (48) 

We will not consider cases for which this equation has no roots since it can be 
shown that the field equations are then not consistent. As a result according to 
the field equations, the theory reduces to General Relativity with a cosmological 
constant given by: 

A = ^ (49) 

For more details see [14] . Therefore investigation of gravitational waves in this 
case is not different from the standard gravity and thus the gravitational wave 
velocity is just cq. 

Now we are going to calculate the gravitational wave velocity for Ricci 
squared gravity. For Ricci squared theory, again we assume that we have not 
negative powers of Ricci squared. That is: 

f(n ab n ab ) = ai n ab n ab + a 2 (n ab n ab ) 2 + ■■■ (50) 

In order to linearize the equations (fT3"]) and (114[) we begin by assuming that the 
metric differs only slightly from the Minkowskian metric, that is 

9ab = Vab + e9ab + 0{e 2 ) (51) 

where e is a small dimensionless parameter and, throughout, we shall neglect 
terms of second order or higher in e. One can easily verify that in the zeroth 
order the affine connection is zero, and thus 

T a bc = + ef a bc + O(e 2 ) (52) 

Neglecting terms of 0(e 2 ), the Ricci tensor is 

K ab = e(f c abtC - T c aC}b ) (53) 

and the Ricci scalar is 

n = g cd n cd = e^ d {t f cdJ -f{ fd ) (54) 

Finally the linearized modified Einstein equation tfTB"]) gives (up to 0(e 2 ) correc- 
tions): 

Kb, c - K c m - lvabv cd (r f cdJ - H f , d ) = o (55) 
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We may seek for wave solutions of the above equation by assuming the 
relations 

r a bc = C bc e ik -* (56) 

h ab = Qab e ik ' x (57) 

where !; bc is symmetric in lower indices because of the fact that we have assumed 
space-times without torsion. Putting equation (|56[) into equation (|55p . we get 
the following gauge condition: 

kcCb ~ hCc = (58) 
The wave equation is derived via linearization of equation Q14p . It gives: 

de(V=99 ab ) + et a ec V cb + sf b ec r, ac + 2a 1 ei 1 ac r 1 bd d e {f{ df ~ f{ f d ) = (59) 

Using the equation V e (y/— gg ab ) = 0, we can simplify the first term and thus 
we have: 

v cbfa c + v acf, ^ + 2ai r 1 ac r 1 bd d e {t{ dJ - £{ f4 ) = (60) 
By using the equations (|56"|) and ([57)) and the above gauge condition, one gets: 

ik c c b =v ad e cd +v bd c d (si) 

After multiplying both sides of the above equation by —ik c , we find out that 
the wave equation becomes 

k 2 C b = -ik c (r 1 ad e cd + V bd S a cd ) (62) 

Taking into account the fact that Palatini theory is invariant under general 
covariance and the so-called projection transformation defined as T% c — > T% c + 
S^Ac + <5°Afc with At, an arbitrary four vector, we have an 16-fold freedom. 
Therefore choosing k c ^ c = 0, we get k 2 — 0. As a result the gravitational wave 
velocity in the nonlinear Ricci squared theory is equal to cq. 

4.3 The electromagnetic wave velocity 

Let us now consider the case of electromagnetic waves by choosing: 

1 



A m = --j d*x^g~F ab F ab (63) 

where the electromagnetic field is defined as 

F ab = D a A b - D b A a = V a A b - V b A a = d a A b - d b A a (64) 
The field equations are thus: 

V a F ab = (65) 
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To find the electromagnetic wave velocity, we consider a solution of type A a — 
e a exp(i/c ■ x). In order to distinguish between the gravitational effects from the 
varying speed of light let's to consider the wave velocity in the two possible local 
frames. 

In the first local frame in which the metric is Minkowskian, for both Ricci 
scalar and Ricci squared gravity we get: 

e • k = (66) 

and 

k 2 = (67) 

The first equation represents gauge invariance and it is equivalent to V ' a A a = 
and the second shows that cem = co in this frame. But remember that in this 
frame the gravitational effects are present. 

In the second local frame for which the connection is zero, there is a distinc- 
tion between Ricci scalar and Ricci squared gravity. For Ricci scalar model one 
has: 

e a (k a + 2id a \nf) =0 (68) 

and 

k 2 + 2ik a d a lnf = (69) 

The first equation is gauge invariance (V a A a = 0), and the second shows that 
generally k 2 ^ 0. Since the extra term is imaginary and is proportional to the 
wave vector one has both attenuation of the wave and that cem ^ cq. In fact 
this extra term is not a constant and thus the electromagnetic wave velocity is 
varying. 

For Ricci squared gravity one has 

e a (k a + ll c ca )^0 (70) 

and 

k 2 + *l a ac k c = (71) 

Again the first equation represents the gauge invariance and the second one 
introduces a varying speed for the electromagnetic waves. 

In order to see the physical implications of this last equation, let us evaluate 
it for FRW cosmological model. A look at equation ([29")) shows that for such 
a model 7° c = jd c X + ^d c uj. Thus the only non-zero element is 7° = q(t) 
which is some function of time. Assuming a plane wave moving in x direction 
with k = (k + ia)x and using equation (|7"Tj) one can easily calculate the phase 
velocity as: 

/ c E m \ _ 



1 



(72) 



and the group velocity as: 



phase A/l + ^n 



1 

2 k 2 



CEM ) = , - • 2 "a /2 (73) 

/ group ( i i q 2 | 



CO 



■Ik 2 
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and also the attenuation factor as: 



The same result is valid for Ricci scalar model as it is clear from equation (|69|) . 
The change in the velocity of electromagnetic waves affects the way atoms radi- 
ate in the far galaxies, and therefore there is a contribution in the cosmological 
redshift from the varying speed of electromagnetic waves. 

4.4 The space— time causal structure constant 

Since the two local frames are conformally related for the case of Palatini Ricci 
scalar theories, the local causal structure is identical to the Minkowskian case. 
This means that in both frames we cst = co- 
in the case of Palatini Ricci squared, it is clear that in the first local frame 
in which the metric is Minkowskian we have c$t = cq. On the other hand 
in the second local frame the space-time structure constant cst has a varying 
character, because of the fact that metric is not locally Minkowskian. In order 
to see how it looks like, let us to investigate it for a cosmological model. We 
consider the FRW model. For such a simple model, we can follow the relations 
easily. In the Robertson- Walker space-time, we have 

dr 2 

ds 2 = g ab dx a dx b = c 2 dt 2 - a 2 (t)[- — j + r 2 {d9 2 + sin 2 9dcj) 2 )} (75) 

and u a = (1, 0, 0, 0). By using of the equations (I57| and (ITS"]) we can find Q a b 

Goo = A.goo + A(cj - 1) = Xuj 
Gij = Agy ► Qu = Xgu 
Got = Agoi = 

So we have 

Gab = diag (V ~*f® , -Xr 2 a 2 (t), -\r 2 sin 2 9a 2 (t)\ (76) 

In order to go to the first local frame for which g a j, is locally r) a b, one can use 
the relation 

g ab = e1{x)e){x)rf* (77) 
with the local tetrad fields given by: 



I Vl - kr 2 1 1 \ 

&i {x) = diag ^ ___ ___ J (78) 

It is quiet obvious that in this frame: 

cst = cq (79) 
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On the other hand, in the second local frame which is local inertial frame, 
we shall have: 

g ab = £?e^ (80) 
Investigation of equation ([76]) shows that the tetrad fields are given by: 



„ a / 1 VI - kr 2 1 1 \ 

e° = diag ^^=, — — , -p (81) 

V \/\uj VAa(t) VAra(i) \J\rsin6a(t) I 

In order to derive cst in this frame, we have to obtain the form of the space-time 
metric using the relation: 

gabete) = 9lJ (82) 

Therefore: 

g ij = d iag(l-,-l,-l,-l) (83) 
The line element in this frame is thus: 

ds 2 = \[-c 2 dt 2 - dr 2 - r 2 d8 2 - r 2 sin 2 8d<j) 2 } (84) 

Assuming a radial null ray, we can find the signal speed which is cst- 

1 . , 

c S T = —j=CQ (85) 
V w 

For a simple model for which f(7Z a blZ ab ) = ai7Z a b'R- ab we obtain 



/ l + 2oiA 

CST = V TT^ C0 (86) 

Let us consider this equation for three situations; namely for dust, radiation 
and cosmological constant. Using the relation lZ ab lZ a b = A 2 + 3S 2 and the 
equations (f3"4"|) , ([33]) we get: 

A + 3aiA 2 - 3H - 3aiS 2 = 2/tp (87) 

and 

- A - aiA 2 - 3 + aiS 2 = -2np (88) 
For dust one has p = and p ~ a -3 , so 

cst ^ (1 + aiKa~ 3 )c (89) 

And for radiation we have p ~ a -4 and p ~ 5 a ~ 4 i so that 

c S T ^ (1 + ^aiKa~ 4 )c (90) 

Finally for cosmological constant cj will be a constant so that the velocity of 
signal is a constant less that Cq. 

Since the region that particles can have interaction with eachother is defined 
using cst, a varying cst can have observable results on problems like the horizon 
problem. 
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5 Conclusions 



In order to conclude let us to stress again that five different velocities can 
be distinguished when one wants to let the velocity of light to vary. In the 
standard theory all these five velocities are equal to the constant Co = 3 x 
10 s meters/ second. 

The first one appears in the coupling constant of gravity and matter and 
here we chose it to be constant. Calling it ce we have ce = cq. 

The second velocity is the gravitational wave velocity cqw- We assumed 
that no negative powers of the Ricci scalar or Ricci squared is present and then 
we have shown that cgw in both scalar and Ricci squared theories is equal to 
Co- 

Since in the theories we concern here the connection is not the metric con- 
nection, two different local frames are distinguishable. The first one is the one 
in which the metric is locally Minkowskian. We saw that the trajectory of a 
test particle in this frame is a straight line. But this does not simply mean that 
this frame is the local inertial (no-gravity) frame. Because gravity arises as 
coupling of the connection to other physical equations and also via the concept 
of parallelism and in this frame the connection is not zero. In the second local 
frame in which the connection is zero, the test particle trajectory is not a line. 
We saw that one can assign the extra terms in the particle's trajectory equation 
as the force due to variation of the clock synchronization speed, the speed of 
the signal used to make faraway clocks synchronized. 

The fourth speed is the electromagnetic wave velocity. We saw that for 
both Palatini non-linear Ricci scalar and Ricci squared cases, the velocity of 
electromagnetic waves is simply cq in the first local frame. But in this frame 
gravity is present and this velocity is the velocity of electromagnetic waves 
modulated with gravitational effects. On the other hand in the second local 
frame which is the local inertial frame, the velocity of electromagnetic waves is 
not Co- It is varying. 

The fifth velocity appears as the space-time causal structure constant or the 
information velocity. In the Palatini f(lZ) theory since the two local frames 
are conformally related, we have cst = cq in both frames. But for the case of 
Palatini Ricci squared, only in the first local frame we have cst — cq,. In the 
second local frame cst is varying. We have investigated it for FRW cosmological 
model and shown that in this local frame that the information velocity is -^Cq- 

To sum up table |T]) may help. 



At this end it is clear that if we use the metric formulation, cem and cst 
would be constant and equal to cq. This is because of the fact that in the metric 
formulation the connection and the Christoffel symbols are equal and thus there 
is only one local frame. In this local frame the metric is Minkowskian and thus 
c st — c o a nd also Maxwell equations have their special relativistic form so we 
have cem = c Q . 
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Palatini f(lZ) gravity 


Palatini f(n ab TZ a b) gravity 


CE 


co 


CO 


CGW 


CO 


CO 


cem 


varying 


varying 


CST 


c 


varying (^c ) 


cc 


varying (-^=) 


varying 



Table 1: Comparison of different speeds of light in the Palatini f(lZ) and 
f(K ab K ab ) theories. 
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